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Abstract 

We consider the applications of a numerical-analytical approach based on multiscale vari- 
ational wavelet technique to the systems with collective type behaviour described by 
some forms of Vlasov-Poisson/Maxwell equations. We calculate the exact fast convergent 
representations for solutions in high-localized wavelet-like bases functions, which corre- 
spond to underlying hidden (coherent) nonlinear eigenmodes. This helps to control stabil- 
ity/unstability scenario of evolution in parameter space on pure algebraical level. 
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Abstract 



or self-similar decompostion [2] like 



We consider the applications of a numerical-analytical 
approach based on multiscale variational wavelet tech- 
nique to the systems with collective type behaviour de- 
scribed by some forms of Vlasov-Poisson/Maxwell equa- 
tions. We calculate the exact fast convergent represen- 
tations for solutions in high-localized wavelet-like bases 
functions, which correspond to underlying hidden (coher- 
ent) nonlinear eigenmodes. This helps to control stabil- 
ity/unstability scenario of evolution in parameter space on 
pure algebraical level. 



1 INTRODUCTION 

In this paper we consider the applications of numerical- 
analytical approach based on multiscale variational wavelet 
technique to the systems with collective type behaviour de- 
scribed by some forms of Vlasov-Poisson/Maxwell equa- 
tions [1], [2]. Such approach may be useful in all mod- 
els in which it is possible and reasonable to reduce all 
complicated problems related with statistical distributions 
to the problems described by the systems of nonlinear or- 
dinary/partial differential/integral equations with or with- 
out some (functional) constraints. In periodic accelera- 
tors and transport systems at the high beam currents and 
charge densities the effects of the intense self-fields, which 
are produced by the beam space charge and currents, de- 
terminine (possible) equilibrium states, stability and trans- 
port properties according to underlying nonlinear dynamics 
[2]. The dynamics of such space-charge dominated high 
brightness beam systems can provide the understanding of 
the instability phenomena such as emittance growth, mis- 
match, halo formation related to the complicated behaviour 
of underlying hidden nonlinear modes outside of perturba- 
tive tori-like KAM regions. Our analysis is based on the 
variational-wavelet approach from [3]-[17], which allows 
us to consider polynomial and rational type of nonlinear- 
ities. In some sense our approach is direct generaliztion 
of traditional nonlinear 6F approach in which weighted 
Klimontovich representation 
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where s(x — Xji) is a shape function of distributing par- 
ticles on the grids in configuration space, are replaced by 
powerful technique from local nonlinear harmonic analy- 
sis, based on underlying symmetries of functional space 
such as affine or more general. The solution has the mul- 
tiscale/multiresolution decomposition via nonlinear high- 
localized eigenmodes, which corresponds to the full mul- 
tiresolution expansion in all underlying time/phase space 
scales. Starting from Vlasov-Poisson equations in part 2, 
we consider the approach based on multiscale variational- 
wavelet formulation in part 3. We give the explicit repre- 
sentation for all dynamical variables in the base of com- 
pactly supported wavelets or nonlinear eigenmodes. Our 
solutions are parametrized by solutions of a number of re- 
duced algebraical problems one from which is nonlinear 
with the same degree of nonlinearity as initial problem and 
the others are the linear problems which correspond to the 
particular method of calculations inside concrete wavelet 
scheme. Because our approach started from variational for- 
mulation we can control evolution of instability on the pure 
algebraical level of reduced algebraical system of equa- 
tions. This helps to control stability/unstability scenario 
of evolution in parameter space on pure algebraical level. 
In all these models numerical modeling demonstrates the 
appearance of coherent high-localized structures and as a 
result the stable patterns formation or unstable chaotic be- 
haviour. 

2 VLASOV-POISSON EQUATIONS 

Analysis based on the non-linear Vlasov equations leads 
to more clear understanding of collective effects and non- 
linear beam dynamics of high intensity beam propagation 
in periodic-focusing and uniform-focusing transport sys- 
tems. We consider the following form of equations (ref. 
[1] for setup and designation): 
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The corresponding Hamiltonian for transverse single-par- 
ticle motion is given by 

H(x,y,p x ,p y ,s) = hpl + p 2 y ) + -[k x (s)x 2 (6) 

+k y (s)y 2 } + Hi(x,y,p x ,p y ,s) + ip(x,y,s), 

where Hi is nonlinear (polynomial/rational) part of the 
full Hamiltonian and corresponding characteristic equa- 
tions are: 

— + k x (s)x + — ip(x,y, s) = (7) 
( ^ + k y {s)y + -^iP(x,y,s) = (8) 

3 MULTISCALE REPRESENTATIONS 

We obtain our multiscale/multiresolution representations 
for solutions of equations (3)-(8) via variational-wavelet 
approach. We decompose the solutions as 

oo 

fb(s,x,y,p x ,p y ) = ^2 ®3 l f( s , x ,y,Px,Py) (9) 

i—i c 

oo 

ip(s, x,y) = Y^ ®S j ip(s, x, y) (10) 



x{s) = Y, ®S k x(s), y(s) = ® st y( s ) (11) 

where set (i c ,j c , k c ,£ c ) corresponds to the coarsest level 
of resolution c in the full multiresolution decomposition 
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Introducing detail space Wj as the orthonormal comple- 
ment of Vj with respect to V^+i : Vj+i = Vj (J) Wj, we 
have for /, ip,x,yC i 2 (R) from (9)-(l 1): 



(13) 



In some sense (9)-(ll) is some generalization of the old 
5F approach [1], [2]. Let L be an arbitrary (non) line- 
ar differential/integral operator with matrix dimension d, 
which acts on some set of functions ^> = ^(s,x) = 

(V(s, x), . . . , f d (s, x)),s,xenc R n+1 from L 2 (n): 



= L(R(s,x),s,x)V(s,x) = 0, (14) 

(x are the generalized space coordinates or phase space co- 
ordinates, s is "time" coordinate). After some anzatzes [3]- 
[17] the main reduced problem may be formulated as the 
system of ordinary differential equations 

Qi(f)¥ i = Pi(f,s), / = (A, •",/»), (15) 
ds 

i = 1, ... ,7i, maxcfeg Pi — p, max deg Qi = q 



or a set of such systems corresponding to each indepen- 
dent coordinate in phase space. They have the fixed initial 
(or boundary) conditions fi(0), where Pi, Qi are not more 
than polynomial functions of dynamical variables fj and 
have arbitrary dependence on time. As result we have the 
following reduced algebraical system of equations on the 
set of unknown coefficients X k of localized eigenmode ex- 
pansion (formula (17) below): 



L{Q ih \,ai) = M(P ij ,\,pj), 



(16) 



where operators L and M are algebraization of RHS and 
LHS of initial problem ( p"5| ) and A are unknowns of re- 
duced system of algebraical equations (RSAE) (|T6|). Af- 
ter solution of RSAE (|T(]) we determine the coefficients 
of wavelet expansion and therefore obtain the solution of 
our initial problem. It should be noted that if we consider 
only truncated expansion with N terms then we have from 
( [To] ) the system of N xn algebraical equations with degree 
£ — max{p, q} and the degree of this algebraical system 
coincides with degree of initial differential system. So, we 
have the solution of the initial nonlinear (rational) problem 
in the form 



N 



(17) 



where coefficients A^" are the roots of the corresponding re- 
duced algebraical (polynomial) problem RSAE ( 16 ). Con- 
sequently, we have a parametrization of solution of ini- 
tial problem by the solution of reduced algebraical prob- 
lem (JTo]). The obtained solutions are given in the form 
(|l7|), where ft (t) are basis functions obtained via multires- 
olution expansions (9)-(ll), (13) and represented by some 
compactly supported wavelets. As a result the solution of 
equations (3)-(8) has the following multiscale/multireso- 
lution decomposition via nonlinear high-localized eigen- 
modes, which corresponds to the full multiresolution ex- 
pansion in all underlying scales (13) starting from coars- 
est one (polynomial tensor bases are introduced in [17]; 
x = (x,y,p x ,p y )): 



*(s,x)= Y ay' 11 ' ® V j (s,x), 



(18) 
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Formula ( p"8| ) gives us expansion into the slow part ^ 
and fast oscillating parts for arbitrary N, M. So, we may 
move from coarse scales of resolution to the finest one 
for obtaining more detailed information about our dynam- 
ical process. The first terms in the RHS of formulae (18) 
correspond on the global level of function space decom- 
position to resolution space and the second ones to detail 
space. It should be noted that such representations give 
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Figure 1 : Eigenmode of level 1 . 




Figure 2: Stable waveleton pattern. 



the best possible localization properties in the correspond- 
ing (phase)space/time coordinates. In contrast with differ- 
ent approaches formulae (18) do not use perturbation tech- 
nique or linearization procedures. So, by using wavelet 
bases with their good (phase) space/time localization prop- 
erties we can describe high-localized (coherent) structures 
in spatially-extended stochastic systems with collective be- 
haviour. Modelling demonstrates the appearance of stable 
patterns formation from high-localized coherent structures 
or chaotic behaviour. On Fig. 1 we present contribution 
to the full expansion from coarsest level (waveleton) of 
decomposition (18). Fig. 2, 3 show the representations 




for full solutions, constructed from the first 6 eigenmodes 
(6 levels in formula (18)), and demonstrate stable local- 
ized pattern formation and chaotic-like behaviour outside 
of KAM region. We can control the type of behaviour on 
the level of reduced algebraical system (16). 
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Figure 3: Chaotic-like behaviour. 
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